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Abstract 

A typical phase field approach for describing phase separation and coarsening phenomena in alloys 
is the Cahn-Hilliard model. This model has been generalized to the so-called Cahn-Larche system by 
combining it with elasticity to capture non-neglecting deformation phenomena, which occur during 
phase separation and coarsening processes in the material. In order to account for damage effects, 
we extend the existing framework of Cahn-Hilliard and Cahn-Larche systems by incorporating an 
internal damage variable of local character. This damage variable allows to model the effect that 
damage of a material point is influenced by its local surrounding. The damage process is described 
by a unidirectional rate-dependent evolution inclusion for the internal variable. For the introduced 
Cahn-Larche systems coupled with rate-dependent damage processes, we establish a suitable notion 
of weak solutions and prove existence of weak solutions. 
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1 Introduction 

Due to the ongoing miniaturization in the area of micro-electronics the demands on strength and lifetime 
of the materials used is considerably rising, while the structural size is continuously being reduced. 
Materials, which enable the functionality of technical products, change the microstructure over time. 
Phase separation and coarsening phenomena take place and the complete failure of electronic devices like 
motherboards or mobile phones often results from micro-cracks in solder joints. 

Solder joints, for instance, are essential components in electronic devices since they form the electrical 
and the mechanical bond between electronic components like micro-chips and the circuit-board. The 
Figures and [^illustrate the typical morphology in the interior of solder materials. At high temperatures, 
one homogeneous phase consisting of different components of the alloy is energetically favourable. If the 
temperature is decreased below a critical value a fine microstructure of two or more phases (different 
compositions of the components of the material) arises on a very short time scale. The formation of 
microstructures, also called phase separation or spinodal decomposition, take place to reduce the bulk 
chemical free energy. Then coarsening phenomena occur, which are mainly driven by decreasing interfacial 
energy. Due to the misfit of the crystal lattices, the different heat expansion coefficients and the different 
elastic moduli of the components, very high mechanical stresses occur preferably at the interfaces of the 
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phases. These stress concentrations initiate the nucleation of micro-cracks, whose propagation can finally 
lead to the failure of the whole electronic device. 





Figure 1: Left: Solder ball and micro-structural coarsening in eutectic Sn-Pb; Right: a) directly after 
solidification, b) after 3 hours, and c) after 300 hours |HCW91j : 



Figure 2: Initiation and propagation of cracks along the phase boundary |FBFD06] . 

The knowledge of the mechanisms inducing phase separation, coarsening and damage phenomena is of 
great importance for technological applications. A uniform distribution of the original materials is aimed 
to guarantee evenly distributed material properties of the sample. For instance, mechanical properties, 
such as the strength and the stability of the material, depend on how finely regions of the original 
materials are mixed. The control of the evolution of the microstructure and therefore of the lifetime of 
materials relies on the ability to understand phase separation, coarsening and damage processes. This 
shows the importance of developing reliable mathematical models to describe such effects. 

In the mathematical literature, coarsening and damage processes are treated in general separately. 
Phase separation and coarsening phenomena are usually described by phase-field models of Cahn-Hilliard 
type. The evolution is modeled by a parabolic diffusion equation for the phase fractions. To include elastic 
effects, resulting from stresses caused by different elastic properties of the phases, Cahn-Hilliard systems 
are coupled with an elliptic equation, describing the quasi-static balance of forces. Such coupled Cahn- 
Hilliard systems with elasticity are also called Cahn-Larche systems. Since in general the mobility, stiffness 
and surface tension coefficients depend on the phases (see for instance |BDM07) and |BDDM07] for the 
explicite structure deduced by the embedded atom method), the mathematical analysis of the coupled 
problem is very complex. Existence results were derived for special cases in ICarOni ICMWinl iBPOb] 
(constant mobility, stiffness and surface tension coefficients), in [BCD^O^ (concentration dependent 
mobility, two space dimensions) and in [PZ08) in an abstract measure-valued setting (concentration 
dependent mobility and surface tension tensors). For numerical results and simulations we refer [WeiOIl 
IMctOSIIBMTO] . 

Damage models for elastic materials have been analytically investigated for the last ten years. In 
the simplest case, the damage variable is a scalar function and describes the local accumulation of 
damage in the body. The damage process is typically modeled as a unidirectional evolution, which 
means that damage can increase, but not decrease. Based on the model developed in |FN96| . the damage 
evolution is described by an equation of balance for forces which is coupled with a unidirectional parabolic 
[BSS051IFK091 IGia05| or rate-independent |MR06l IMRZlOj evolution inclusion for the damage variable. 
The models studied in [FK091 IMROBl IGia05) also include the effect that the applied forces have to pass 
over a threshold before the damage starts to increase. 
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In this work, we introduce a mathematical model describing both phenomena, phase separation/coars¬ 
ening and damage processes, in a unifying model. We focus on the analytical modeling on the meso- and 
macroscale. To this end, we couple phase-field models of Cahn-Larche type with damage models. The 
evolution system consists of an equation of balance for forces which is coupled with a parabolic evolution 
equation for the phase fractions and a unidirectional evolution inclusion for the damage variable. The 
evolution inclusion also comprises the phenomenon that a threshold for the loads has to be passed before 
the damage process increases. 

The main aim of the present work is to show existence of weak solutions of the introduced model 
for rate-dependent damage processes. A crucial step has been to establish a suitable notion of weak 
solutions. We first study the model with regularization terms and prove existence of weak solutions for 
the regularized model based on a time-incremental minimization problem with constraints due to the 
unidirectionality of the damage. The regularization allows us to prove an energy inequality which occurs 
in the weak notion of our coupled system. The major task has been to prove convergence of the time 
incremental solutions for the regularized model when the discretization fineness tends to zero. In this 
context, several approximation results have been established to handle the damage evolution inclusion 
and the unidirectionality of damage processes. More precisely, the internal variable z, describing damage 
effects, is bounded with values in [0,1] and monotonically decreasing with respect to the time variable. 
The main results are stated in Sections |4.1| and |4.2[ see Theorems |4.4| and |4.6[ 

To the best of our knowledge, phase separation processes coupled with damage are not studied yet 
in the mathematical literature. However, promising simulations were carried out in the context of phase 
field models of Cahn-Hilliard and Cahn-Larche type with damage, see |USG071 lGIJaMM+07 . 

The paper is organized as follows: We start with introducing a phase field model of Cahn-Larche type 
coupled with damage, cf. Section Then we state some assumptions for this model, see Section In 
Section we establish a suitable notation for weak formulations of solutions for the introduced model 
and a regularized version of the model and state the main results. Section [5.2| is devoted to the existence 
proof for the regularized Cahn-Larche system coupled with damage. Finally, we pass to the limit in the 
regularized version, which shows the existence of weak solutions of the original model, see Section |5.3| 

2 Model 

We consider a material of two components occupying a bounded Lipschitz domain H C . The state of 
the system at a fixed time point is specified by a triple q = {u, c, z). The displacement field rt : H 
determines the current position x -I- u(x) of an undeformed material point x. Throughout this paper, 
we will work with the linearized strain tensor e(u) = ^(Vu -I- (Vu)^), which is an adequate assumption 
only when small strains occur in the material. However, this assumption is justified for phase-separation 
processes in alloys since the deformation usually has a small gradient. The function c : H K is a phase 
field variable describing a scaled concentration difference of the two components. To account for damage 
effects, we choose an isotropic damage variable z : H —>■ K, which models the reduction of the effective 
volume of the material due to void nucleation, growth, and coalescence. The damage process is modeled 
unidirectional, i.e. damage may only increase. Self-healing processes in the material are forbidden. No 
damage at a material point x S H is described by z{x) = 1, whereas z{x) = 0 stands for a completely 
damaged material point x S H. We require that even a damaged material can store a small amount of 
elastic energy. Plastic effects are not considered in our model. 

2.1 Energies and evolutionary equations 

Here, we qualify our model formally and postpone a rigorous treatment to Section The presented 
model is based on two functionals, i.e. a generalized Ginzburg-Landau free energy functional £ and a 
damage dissipation potential TZ. The free energy density ip of the system is given by 

r 

(^(e,c,Vc,z,Vz) := ^|Vc|2 + -|Vzr + Wch(c) + Wei(e,c,z), 7, > 0, (1) 
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where the gradient terms penalize spatial changes of the variables c and z, W^h denotes the chemical 
energy density and Wei is the elastically stored energy density accounting for elastic deformations and 
damage effects. For simplicity of notation, we set 7 = <5 = 1. 

The chemical free energy density Wch has usually the form of a double well potential for a two phase 
system. For a rigorous treatment, we need the assumptions (A1)-(A 61 , see Section]^ Hence, in particular, 
classical ansatzes such as 

Wch = (1 - C2)2 


fit in our framework. 

The elastically stored energy density Wei due to stresses and strains, which occur in the material, is 
typically of quadratic form, i.e. 


Wei(c,e) = ^(e-e*(c)) : C(c)(e - e*(c)). 


( 2 ) 


Here, e*(c) denotes the eigenstrain, which is usually linear in c, and C(c) € is a fourth order 

stiffness tensor, which is symmetric and positive definite. If the stiffness tensor does not depend on the 
concentration, i.e. C(c) = C, we refer to homogeneous elasticity. 

To incorporate the effect of damage on the elastic response of the material, Wei is replaced by 

Wei = (^(^) + V) Wei, (3) 


where $ : [ 0 , 1 ] —>■ M+ is a continuous and monotonically increasing function with d)( 0 ) = 0 and 77 > 0 is 
a small value. The small value 77 > 0 in ([^ is introduced for analytical reasons, see for instance 

Rigorous results in the present work are obtained under certain growth conditions for the elastic 
energy density Wei, see Section]^ These conditions are, for instance, satisfied for Wei as in (§ in the 
case of homogeneous elasticity. 

The overall free energy £ of Ginzburg-Landau type has the following structure: 


£{u, c, z) 


£(u,c,z)+ / /[o,oo) (2:) dx, 


£{u,c,z) 


/ (/ 7 (e(M), c, Vc, z, Vz) dx. 

Jn 


(4) 


Here, d[Q ^ 0 ) signifies the indicator function of the subset [0, 00 ) C M, i.e. d[o,oo)(3:) = 0 for a: G [0, 00 ) and 
-^[ 0 , 00 ) ( 2 :) = 00 for X < 0. We assume that the energy dissipation for the damage process is triggered by 
a dissipation potential TZ of the form 


7 ^(z) :=TZ{z)+ [ /(-oo.o]{2)dx, 

Jn 

'R{z) := f —az + dx for a > 0 and /3 > 0. 

Jn 2 


(5) 


Due to /3 > 0, the dissipation potential is referred to as rate-dependent. In the case /3 = 0, which is 
not considered in this work, TZ is called rate-independent. We refer for rate-independent processes to 
[EM061IMT991IMR061 IMR.ZIOI IRonlO) and in particular to [Mie05) for a survey. 

The governing evolutionary equations for a system state q = {u, c, z) can be expressed by virtue of the 
functionals Q and ([^. The evolution is driven by the following elliptic-parabolic system of differential 
equations and differential inclusion: 


Diffusion : 

dtc= Afi{u,c,z), 

(6a) 

Mechanical equilibrium : 

div[a(e{u),c, z)) = 0, 

(6b) 

Damage evolution : 

0 G dz£{u, c, z) -1- diTZ{dtz), 

(6c) 


where a = a{e, c, z) := 9e<p(e, c, Vc, z, Vz) denotes the Cauchy stress tensor and 77 is the chemical poten¬ 
tial given by 77 = 77(77,0, z) := i9c<p(e, c, Vc, z, Vz) — div(i 9 vc< 75 (e, c, Vc, z, Vz)). Equation ( 6 a) is a fourth 
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order quasi-linear parabolic equation of Cahn-Hilliard type and describes phase separation processes for 
the concentration c while the elliptic equation (6b) constitutes a quasi-static equilibrium for u. This 
means physically that we neglect kinetic energies and instead assume that mechanical equilibrium is at¬ 
tained at any time. The doubly nonlinear differential inclusion (6c) specifies the flow rule of the damage 
profile according to the constraints 0 < z < 1 and dtZ < 0 (in space and time). The inclusion (6c) has to 
be read in terms of generalized sub-differentials. 

We choose Dirichlet conditions for the displacements u on a subset T of the boundary dfl with 
•^n-i(r) > 6 : [0,T] X r — K" be a function which prescribes the displacements on T for a fixed 

chosen time interval [0,T]. The imposed boundary and initial conditions and constraints are as follows: 


Boundary displacements : u(t) = b(t) on T for all t S [0,r], (IBCf) 

Initial concentration: c(0) = c^ in (IBC2) 

Initial damage : 0 < 2:(0) = < 1 in 17, (IBC3) 

Damage constraints : 0 < z < 1 and dfZ < 0 in Sir. (IBC4) 


Moreover, we use homogeneous Neumann boundary conditions for the remaining variables on (parts of) 
the boundary: 


cr • 1 / = 0 on 917 \ r. 

(IBC5) 

V^(t) • 1 / = 0 on 917, 

(IBC6) 

Vc(t) • 1 / = 0 on 917, 

(IBC7) 

Vz(t) ■ V — 0 on 917, 

(IBC8) 


where v stands for the outer unit normal to 917. 

We like to mention that mass conservation of the system follows from the diffusion equation (6a) and 
(p^, i.e. 


j c{t) — c° dx = 0 for all t G [0, T], 
Jn 


3 Assumptions and Notation 

In the following, we collect all assumptions and constants which are used for a rigorous analysis in the 
subsequent sections. 

(i) Setting. 17 C M” is a bounded domain with Lipschitz boundary, n G {1,2,3}, p > n, (3 > 0, 

Wei e X M X K;K+), Wch G C'i(K;M+), Wei(e,c,z) = Wei(e*,c,z) for all e G and 

c, z G M. Furthermore, C > 0 always denotes a constant, which may vary from estimate to estimate, 
and [0,T] is the time interval of interest. 

(ii) Convexity and growth assumptions. The function Wei is assumed to satisfy for some constants 


T] > 0 and C > 0 the following estimates: 

? 7 |ei - e 2 p < (9eWei(ei,c,z) - 9 eWei(e 2 , c, z)) : (d - 62 ), (Al) 

W;i(e,c,z)<C(|ep + |c|2 + l), (A2) 

|9eWei(ei -h e 2 ,c, z)| < C'(Wei(ei,c, z) -t |e 2 | -b 1), (A3) 

|9eWei(e,c,z)| <C(|e| + |cp + l), (A4) 

\d^We\{e, c, z)| < Cdep -b |cp -b 1) (A5) 

for arbitrary c G M, z G [0,1] and symmetric e, ei, 62 G K”^". 

The chemical energy density function Wch satisfies 

|9cWch(c)| <C(|c|2*/2 + i) (A6 ) 
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for sonie constant C > 0. For dimension n = 3, the constant 2* denotes the Sobolev critical 

In the two dimensional case n = 2, the constant 2* can be an arbitrary 


exponent given by . 


positive real number and in one space dimension (A 6 ) can be dropped. 


(hi) Boundary displacements. We assume that F is a ^-measurable subset of dil with 'H”“^(F) > 0 
and that the boundary displacement 6 : [0, T] x F —>■ K” may be extended by & S T]; W^’°° 

(17; K”)) such that &(t)|r = Ki)|r in the sense of traces for a.e. t G [0, T]. In the following, we write 
b instead of b. 


Remark 3.1 Conditions (All, 


and (A31 imply the following estimates 


|aetFei(e,c,z)| <C(|e| + |c|2 + I), 

rylep - C{\c\^ + 1) < Wei(e, c, z) 

for some appropriate constants rj > 0 and C > 0, cf. ^GarOfA Section 3.2] for dllbt . 


(I la) 
( 11 b) 


We introduce some auxiliary spaces to shorten the notation for the construction of solution 
the evolutionary problem. First of all, we define the trajectory space Q for the limit problem 
as 


curves of 


r mg L°°([0,T];Ri(f7;K")), 

Q:= lq={u,c,z) with c G L°° {[0,T]-, H\n)) n H\[0,T], 

[ z G L'^{[0,T];W^'P{n)) n H^i[0,T]-,L^{n)) 


Based on Q, the set of admissible functions of the viscous problem (see Section]^ is 

Q'' ■= {q={u,c,z) G Q\cGH\[0,T]-,L^{n)) and m G L°°([0,T]; 1F1’^(17;K”))}. 


It will be convenient for the variational formulation to define Sobolev spaces with functions taking 
only non-negative and non-positive values, respectively, and Sobolev spaces consisting of functions with 
vanishing traces on the boundary F: 

:= {C G lF^’’'(f7) I C > 0 a.e. in 17}, 

Wy(17) := {C G 1F^’’'(17) I C < 0 a.e. in 17), 
lFp’'^(17;K"’) ;= G 1F^’’’(17; M") | = 0 in the sense of traces} 

for r G [l,oo]. In this context, : IF^’’'(17) —>• {0, oo} denote the indicator functions given by 





if C G 1A1’"(17), 
else. 


Since Cahn-Hilliard systems can be expressed as H ^-gradient flows, we introduce the following spaces 
in order to apply the direct method in the time-discrete version (see Section]^: 

Fo := |c G R^(17) I J Cda; = 0 

Fo:= {C G (Ri(17))* I (C,l)(^,).,^, =0}. 

This permits us to define the operator (—A)“^ : Vq ^ Vq as the inverse of the operator —A : Vq ^ Vq, u i-G 
(Vu, V-)j;, 2 (Q). The space Fq will be endowed with the scalar product {u,v)y^ := (V(—A)“^u, V(—A)“^M)i 2 (f 2 ). 

We end this section by introducing some notation which is frequently used for some approximation 
features in this paper. The expression Bji{K) denotes the open neighborhood with width i? > 0 of a 
subset K C K". Whenever we consider the zero set of a function C G IF^’^(17) for p > n abbreviated in 
the following by {C = 0} we mean {x G 17 | (]{x) = 0} by taking the embedding IF^’*’(17) ^ C^(fl) into 
account. We adapt the convention that for two given functions Ci? G L^([0, T]; IF^’*’(17)) the inclusion 
{C = 0} 2 {C = 0} is an abbreviation for {C(t) = 0} A {^(t) = 0} for a.e. t G [0,r]. 
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4 Weak formulation and existence theorems 


Existence results for multi-phase Cahn-Larche systems without considering damage phase fields are shown 
in [GarOO] provided that the chemical energy density Wch can be decomposed into -I- with convex 
and linear growth behavior of dcW^^ (see |Gar001 Section 3.2] for a detailed explanation). Logarithmic 
free energies Wch are also studied in |GarOO| as well as in |Gar05bj . Further variants of Gahn-Larche 
systems are investigated in jCMPOO) . [BP05j . |BCD+02] and |Gar05aj . 

Purely mechanical systems with rate-independent damage processes are analytically considered and 
reviewed for instance in [MR06| and [MRZlOj . The rate-independence enables the concept of the so-called 
global energetic solutions (see Remark 4.2 (i)) to such systems. 


Goupling rate-independent systems with other (rate-dependent) processes (such as with inertial or 
thermal effects) may lead, however, to serious mathematical difficulties as pointed out in [RoulOj . 

In our situation where the Gahn-Larche system is coupled with rate-dependent damage, we will treat 
our model problem analytically by a regularization method that gives better regularity property for c 
and integrability for u in the first instance. A passage to the limit will finally give us solutions to 
the original problem. In doing so, the notion of a weak solution consists of variational equalities and 
inequalities as well as an energy estimate, inspired by the concept of energetic solutions in the framework 
of rate-independent systems. 


4.1 Regularization 

The regularization, we want to consider here, is achieved by adding the term eAdtC to the Cahn-Hilliard 
equation (referred to as viscous Gahn-Hilliard equation |BP05j l and the 4-Laplacian ediv(|VMpVu) to 
the quasi-static equilibrium equation of the model problem. The classical formulation of the regularized 
problem for e > 0 now reads as 


dtc = A{-Ac + dcWch{c) + dcWeiie{u),c,z) +edtc), 
div(tT(e(u), c, z)) + ediv(|VupVw) = 0, 

0 e d^£e{u, c, z) + diTZ{dtz) 

with the regularized energies 


(12a) 

(12b) 

(12c) 


£s{u, c, z) := £{u,c, z) + e / -|VM|^da;, 

Jn 4 


£e{u,c, z) := £{u,c, z) + e / --\^u\ dx. 

Jn 4 


In the following, we motivate a formulation of weak solutions of the system (12al-(12b) admissible for 


curves q = (u,c,z) € Q''. For every t S [0,T], equation (12a) can be translated with the boundary 
conditions in a weak formulation as follows: 


f {dtc{t))(^ dx = — [ V/r(t)-VCd 2 : 

Jn Jn 


(13) 


for all C, € and 


/ /j,(t)Cdx= / Vc(t)-VC+ dcWch(c(t))C+ dcWei(e(u(t)),c(t),z(t))C + s(dtc(t))Cdx (14) 
/n Jn 


for all () S In the same spirit, we rewrite (12b I as 


deWei{e{u{t)),c{t),z{t)) : e(C) -b e\Wu{t)\‘^Wu{t) : V() dx = 0 


( 15 ) 


for all C G lTp’^(r2 ;K") by using the symmetry condition 


5eWei(e, c, z) = {dcWciie, c, z)f for e G 


sym 


, c, z G 
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following from the assumptions in 


Section]^ (i). The differential inclusion (12c) is equivalent to 


0 = d^Se{u{t), c{t),z(t)) + r{t) + di'R.{dtz{t)) + s{t) 

with some r{t) S dly^i.v and s[t) G dl^^i,p^^^{dtz{t)) (see (|^ and ([^ for the definitions of E and 
1Z). This can be expressed to the following system of variational inequalities: 


{dtz{t)) 


+ r{t) + din{dtz{t))X - dtz(t)^ < 


for C G W^’P{n), 
for C G W^’P{n). 


Here, (•, •) denotes the dual pairing between (lT^’^(f2))* and IT^’^(H). This system is, in turn, equivalent 
to the inequality system 


z{t) > 0 and dtz{t) < 0 , 


(16a) 

d^4(9(0) + + din{dtz{t)),dtz{t)'^ > 0 , 


(16b) 

(d:,£e{q{t)) + r{f) + di'R.{dtz{t))x') > 0 

for C G Wb^{n), 

(16c) 

{r{t)X- z{f)) < 0 

for c G wl’P{n). 

(16d) 


Due to the lack of regularity of q, { 16b I cannot be justified rigorously. To overcome this difficulty, we use 
a formal calculation originating from energetic formulations introduced in |MT99| . 


Proposition 4.1 (Energetic characterization) Let q G Q''nC‘^{Llj 


' X K X K) be a smooth solution 


of (13)-(15) with (IBCl )-(|IBC 8 ). Then the following two conditions are equivalent: 


(i) (16bl with r[t) G for all t G [0,T], 

(ii) for all 0 < ti < t 2 < T: 


rt2 






£e{q{t2))+J {diTZ{dtz),dtz) ds 
rt2 

JQ. 


|V/rp + e\dtc\^ dxds - £e{q{ti)) 


< 


deWei{e{u), c, z) : e{dtb)dxds + e 


ri2 


iVupVu : VSt^dxds. 


Proof. We first show for all t G [0,T]: 

{r,dtz{t)) = 0 for all r G dly^i,p^p^^{z{t)). 


(17) 


(18) 


The inequality 0 < {r,dtz{t)) follows directly from (16d) by putting C, = z{t) — dtz(t). The ’>’ - part 
can be shown by an approximation argument. Applying Lemma 5.1 with /m = z{t) and / = z(t) 
and ( = —dtz{t), we obtain a sequence {Cm} C 1 T(J:’^(H) and constants um > 0 such that —Cm —>■ 
dtz{t) in W^’P{Ll) as M —>• oo and 0 < z{t) — vmCm a.e. in fl for all M G N. Testing (16d) with C = 
z{t) — vmCm shows (r, —Cm) < 0 . Passing to M —>• oo gives (r, dtz{t)) < 0. 


To {ii) => (i) : We remark that (141 and (15) can be written in the following form: 

[ T{t)Ci - £{dtc{t))(i da; = (dc4(g(t)), Ci), 

Jn 

{du£e{q{t))X2) = 0 , 


for all t G [0,T], all Ci G H^{El) and all C 2 G M"). 


(19a) 

(19b) 

























Let to S [0,T). It follows 


£Mto + h)) £e{q^ ^ f°^\diTi{dtz),dtz) dt + / 


h 

ptQ-\-h 

- J 

J to * 

Letting h\0 gives 




*0 


/'to-\-h n 

to J Q. 


L 


|V/^P + £|9tcp dxdt 


pto-\-h p 

deWei{e{u),c,z) : e{dtb) dxdt + e+ / \Vu\^Vu : Vdtbdxdt. 

! Jto Jn 


^^£e{q{to)) + {diii{dtz{to)),dtz{to)) + J |V^(to)|^ + e|9tc(to)Pdx 

< [ deWei{e{u{to)),c{to),z{to)) : e{dtb{to))dx + e [ |Vu(to)pVM(to) : V5t6(to) da; 
JQ Ja 

= {du£s{q{to)),dtb{to)). 


Using the chain rule and (13)-(151 yield 


—£e{q{to)) = {du£e{.q{to)), dtu{to)) + {dc£e{q{to)), dtc{to)) +{d:,£e{q{to)),dtz{to)) 

at s. ^ ✓ s, ^ ✓ 


apply l |19b( 

= {du£e{q{to)),dtbito)} 


apply ( |19a[ l and l |13^ 

-|VAi(to)P - e\dtc{to)\'^ dx + {d^£e{q{to)),dtz{to)). 


In consequence, property (i) follows together with (18). The case to = T can be derived similarly 
by considering the difference quotient of to and to — h. 


To {i) => (ii) : This implication follows from the relation £e{q{t 2 )) — £e{q{ti)) = ^£e{q{t)) dt as 

well as the equations (13)-(15) and (18). ■ 


Remark 4.2 


(i) In the rate-independent case /3 = 0 and for convex £e with respect to z, condition (16c) can be 
characterized by a stability condition which reads as 


£e{u{t), c{t),z{t)) < £eiu{t),c{t),C) + 'R-iC “ z{t)) 


( 20 ) 


for all t € [0,T] and all test-functions C. G lU]:’^(fl). Thereby, ( [T7| ) and ( p0| ) give an equivalent 
description of the differential inclusion { 12c) for smooth solutions. This concept of solutions is 
referred to as global energetic solutions and was introduced in JMTQSff . We emphasize that the 
damage variable z in the rate-independent case /3 = 0 is a function of bounded variation and is 
allowed to exhibit jumps. For a comprehensive introduction, we refer to \AFP0(lf . To tackle rate- 
dependent systems and non-convexity of £e with respect to z, we can not use formulation (20) (cf. 
IMRSOfA \ME.Z10I ). 

(ii) For smooth solutions q, satisfying (13)-(15), the energy inequality 0 and the variational inequality 
(16c), we even obtain the following energy balance: 


£e(.q{t2))+ j {di'R.{dtz),dtz) ds F 

Jti 

= £eiqiti))+ [ 

J t\ 

for all 0 < ti < t 2 < T. 


+ e|9tcp da;ds 




deWe\{e{u),c,z) : e(9(6) da;ds + e 


iVupVw : Wdtbdxds 
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This motivates the definition of a solution in the following sense: 


Definition 4.3 (Weak solution - viscous problem) A triple q = {u,c,z) G Q'' with c(0) = 

2 :( 0 ) = , z > 0 and dtz < 0 a.e. in fix is called a weak solution of the viscous system ( 12a)-(|l2^ 


with initial-boundary data and constraints (IBCl )-( IBC8[ ) if it satisfies the following conditions: 

(i) for all C G 

j {dtc)(f dxdt = — I V/i-VCdxdt, (21) 

where p, G T^([0, T]; satisfies for all C. G T^([0, T]; 

( yLC,dxdt= [ \7c-VC + dcWch{c)C-\-dcWei{e{u),c,z)C-\-e{dtc)Cdxdt, (22) 

J 0^1 J Ot" 

(a) for all^G L4([0,T];Wr’^(f];K’^)) 

f i9eld4i(e(M), c, z) : e(C) + e|VupVu : VCdxdt = 0, (23) 

J 

(Hi) for all C G LP{[0,T];Wb^{n)) n L“(flT) 

0< [ \yz\P-^yz-yC + {d,W,i{eiu),c,z)-a + fiidtz))Cdxdt+ [ {r{t),C{t))dt, (24) 

J J 0 

where r G L^{^t) 'G {^,T\,{W^’P{^)Y^ satisfies for all ( G {fl) and for a.e. tG[0,T] 

{r{t),C - z{t)) <0, (25) 


(iv) for a.e. 0 <ti <t 2 <T 

^eiqih)) + [ a{z{ti) - z{t 2 )) dx + f f ^|5tzpda;ds+ / f | + £:|5tcp dxds 

Jq Jti JQ, Jt\ Jq 

< £e{q{ti)) j f deWe\{e{u), c, z) : e(dtb) dxdse f f |VupVu : Vcitfe deeds. (26) 
J t\ J^ Jt\ 'J ^ 

Theorem 4.4 (Existence theorem - viscous problem) Let the assumptions in 5ectzon[^ be satisfied 
and let c° G z^ G W^'P(fil) with 0 < < 1 a.e. in LI and a viscosity factor e G (0,1] be given. 

Then there exists a weak solution q G Of' of the viscous system (12a)-(12c) in the sense of Definition \4.^ 
In addition: 


r = -X{z=o} [dzWe\{e{u), c, z)]’ 


(27) 


where [•]+ is defined by max{0, •}. 


4.2 Limit problem 


Our main aim in this work is to establish an existence result for the system (12a)-(12c) with vanishing 


e-terms, i.e. with e = 0. In the same fashion as in Section 4.1 we introduce a weak notion of (6a|-(6c) as 
follows. 

Definition 4.5 (Weak solution - limit problem) A triple q = {u, c, z) G Q with z(0) = z°, z > 0 
and dtZ < 0 a.e. in LIt is called a weak solution of the system (6a)-(6c) with boundary and initial 
conditions (IBCl )-(IBC8|) if it satisfies the following conditions: 
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(i) for all C € T]; with dtC G L'^i^r) Bind C,[T) = 0 


f {c — c^)dt( dxdt = ( V^-VCdxdi, 

J O'J' 

where /i S ^^([0, T]; satisfies for all ( € ^^([0, T]; H^{D,)) 

/ ^Cda:di= / Vc • VC + i9cW^ch(c)C + 9cVFei(e(zi), c, 2 :)Cdxdi, 

J CItp J 


(ii) for all C, e L^{\f),T];Hl{Sl-W^)) 


[ deWei{e{u),c,z) : e{()dxdt ^ 0, 

J Qt’ 


(Hi) for all C G LP{[0,T]-,Wb^(n)) n 

0 < / \\7z\P~‘^\7z-S/C + d:^Wei{e{u),c,z)( - a( + P{dtz)((dxdt + f {r{t),C{t)) dt, 

J J 0 

where r G satisfies for all ( G W^’^(ri) a77.d /or a.e. t G [0, T] 

{r{t),C- z{t)) < 0, 


(iv) for a.e. 0 <ti <t 2 <T 

£{q{t 2 ))+ f a{z(ti) - z(t 2 )) dx + f f (3\dtz\'^ dxds + f j |V^pda;ds 

t/Q 'J t\ 'J Vl 'J t\ 'J ^ 


< ^( 9 (^ 1 )) 


rt2 


‘ ti JQ 


deWf,i[e{u),c,z) : e{dtb)dxds. 


Theorem 4.6 (Existence theorem - limit problem) Let the assumptions in Section he satisfied 
and let c° G z^ G W^’P{Lt) with 0 < < 1 a.e. in fl be given. Then there exists a weak solution 

q G Q of the system (6a)-([^ in the sense of Definition |/ . 


5 Proof of the existence theorems 


5.1 Preliminaries 


The proof of Theorem 4.4 is based on recursive functional minimization that comes from an implicit Euler 
scheme of the system (12a)-(12c) with respect to the time variable. To obtain from the time-discrete model 


the time-continuous model (12al-(12cl, we need some preliminary results on approximation schemes for 
test-functions, which will be presented in this section. 


Lemma 5.1 (Approximation of test-functions) Let p > n and /, C G with {C = 0} 3 {/ = 

0}. Furthermore, let {/m}mgn C WlJ:’^(f2) be a sequence with fM f in W^’P{Vl) as M ^ 00 . Then, 
there exist a sequence {Cm}mgn C and constants vm > 0, M G N, such that 

(i) Cm C in as M ^ 00 , 

(ii) Cm < C a.e. in fl for all M gN, 

(Hi) vmQm < /m o-.e. in Lt for all M gN. 
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Proof. Without loss of generality we may assume C ^ 0 on 

Let {(5fe} be a sequence with <5^ 0 as /c —)■ oo and (5^ > 0. Define for every k € N the approximation 

function (k € as 


a := [C-4] + , 

where [•]■•■ stands for max{0, •}. Let 0<a<l— ^bea fixed constant. Then (k € due to 

^ C°’“(D). Furthermore, set the constant Rk, fc € N, to 

Rk '■= (^4/||Cllc0"“(O)) ^ 

It follows {(k = 0} D D n S_Rfc({C = 0}) D D n S_Rfc({/ = 0}). without loss of generality we may assume 
n \ Bji^{{f = 0}) 7 ^ 0 for all fc € N. Furthermore, there exists a strictly increasing sequence {Mk} C N 
such that we find for all fc € N: 

/m > rik/2 a.e. on D \ BR^{{f = 0}) for all M > Mk 

with rjk := inf{/(x) | a; € D \ = 0})} > 0, fc € N, (note that /m —> / in as M ^ oo). 

This implies f'fcCfe < /m a.e. on ft for all M > Mk by setting i>k '■= ■'7fe/(2||CI|L“(n)) > 0. The claim 
follows with Cm := 0 and Vk ■= 1 for M € Mi — 1 } and Cm := CSk vm '■= 4 for each 

MG{Mfc,...,Mfc+i-l}, fccN. ■ 


Lemma 5.2 (Approximation of time-dependent test-functions) Let p > n, q > 1 and /, C G 

L'^([0, T]; IF]i’^(D)) with {C = 0} D {/ = 0}. Furthermore, let {/MjMew C L'?([0, T]; IF]i’^(D)) be a 
sequence with fikiit) f{t) in LF^’^(D) as M ^ oo for a.e. t G [0,T]. Then, there exist a sequence 
{CMjMeN C L®([0, T]; IF^’^(D)) and constants VM,t > 0 such that 

(i) Cm in L«([ 0 , T]; LFbP(D)) as M ^ oo, 

(ii) Cm < C for all M gN (in particular {Cm = 0} D {C = 0 }), 

(Hi) VM,tCM{t) < fM{t) a.e. in D for a.e. t G [0,r] and for all M gN. 

If, in addition, C ^ f o.e. in Ht then condition (Hi) can he refined to 
(Hi)’ Cm < fM o.e. in LIt for all M gN. 

Proof. Let {6k} with \ 0 as A: —>■ oo and (5fc > 0 be a sequence and 0<a<l — ^bea fixed constant. 
We construct the approximation functions Cm € ^*([0, T]; VF^’^(D)), M gN, as follows: 

M 

CM(t) :=^Xa^W[C(A)-4] + , (28) 


where Xa^ ■ defined as the characteristic function of the measurable set given by 


4« — 


pk 

pM 


\(u 


M p 
=k+l 


^ if fc < M, 
if fc = M, 


with 

Ptt := {t G [0,T] I n \ Bi^,(,)({/(t) = 0}) 7^ 0 

and fM{t) > Vk{t)/2 a.e. on D \ Bij^p)({/(t) = 0})|, (29) 
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where the functions Rk,Vk '■ K"*" are defined by 


Rk{t) — (^fc/IIC(^)llco.“(n)) ) 

r]k{t) =mi{f{t,x)\x G \ 5^,(4) ({/(t) = 0})}. 


Here, we use the convention Rk{t) := oo for C{t) = 0. Note that 1 < fc < M, are pairwise disjoint 
by construction. 

Consider a t G [0,T] with fnit) fit) in and ({t) ^ 0 with {C(t) = 0} 3 {f{t) = 0}. Let 

iL G N be arbitrary but large enough such that = 0}) ^ 0 holds. It follows the existence 

of an M > K with t G for all M > M. Therefore, for each M > M exists ak> K such that t G 

i.e. CM(t) = iCit) — Thus CM(t) —t C(^) in W^’^i^) as K ^ oo. Lebesgue’s convergence theorem 

shows (i). 

Property (ii) follows immediately from (28). It remains to show (iii). Let M G N be arbitrary. If 


Cuit) = 0 we set i^M,t = 1- Otherwise we find a unique 1 < k < M with t G A\^ and CM(t) = [C(^) ~ 
This, in turn, implies the existence of a > 0 with VM,tC,M < /m (see proof of Lemma 5.1). 

In the case C, < f, we use instead of (29) the set: 


-pk __ 


■■= {tG[0,r]|||/M(0-/(i)llco(n)<4}. 

With a similar argumentation, {Cm} fulfills (i), (ii) and (iii)’. 


Lemma 5.3 Let p > n and f G L^Ap i)(ri;K"), g g L^iil), z G with f ■ Vz > 0 and 

{/ = 0} O {z = 0} a.e.. Furthermore, we assume that 

I f • VC + ffC da; > 0 for all C G Wf^^iLt) with (C = 0} O {z = 0}. 
an 


Then 


//•VC + 5 Cdx> / [ff]+Cdx for allt^GWb^in). 
an a{z=o} 


Proof. We assume z ^ 0 on 0. Let C G Wf'^iLl) be a test-function. For <5 > 0 small enough such that 
O \ Bsi{z = 0}) 7 ^ 0, we define 


Ci := max{C,-z||CIU“C'5 


with the constant 


Cs := inf {z(a;) | a; G O \ H^dz = 0})} > 0. 
We consider the following partition of ft: 

n = Si U Sy U 

with 


El :=n\Hd{z = 0}), 

E| ■.= nnBsi{z = 0}) n (C < -zKh^c^"}, 
E> ■.= nnBsi{z = 0}) n (C > -dlCIU-o.-d- 

By construction, the sequence {C(5}(5 g(o,i] satisfies 

r ( ) = if a: G El U E^, 

\-z(x)||C||l»Q-\ ifa;GE|. 
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In particular, = 0 on {z = 0} for every S S (0,1] and ^,5 ^ in L°°{{z > 0}) as 5 \ 0. By using the 
assumptions, we estimate 


/• VC + ffCda:- 


/{z=0} 


[5]+Cdx 


= [ /•V(C-a)+5(C-a)da:- f [5]+Cda.+ f f-WCs+gCsdx 
J Vl j { 2 = 0 } j n 


>0 


> //•V(C-a)dx+ / 9{.C-Cs)dx 
Jq J { 2 >o} 

= / /•V(C-a)da:+ f J.y{C-Cs)dx+ / /•V(C-a)dx 

l/Sl l/So 


=0 


=0 


+ [ 9{C-Cs)dx 

J{z>0} 

= K\\l-C^" [ f-^zdx+ [ f-\/Cdx +[ 9{C-Cs)dx 
Jsl '—V—^ JsS Jfz>0\ 


>0 


= r < / VC da; 


> 


[ f-VCdx+[ 9{C-C5)dx. 

/e|\{z= 0} /tz>0> 


{z>0} 


The terms on the right hand side converge to 0 as (5 \ 0. 


5.2 Viscous case 


This section is aimed to prove Theorem 4.4 The initial displacement is chosen to be a minimizer 
of the functional u !->■ fe(u,c°,z°) defined on the space -^^ith the constraint M|r = fe( 0 )|r (the 

existence proof is based on direct methods in the calculus of variations - see the proof of Lemma |5.4| 
below). We now apply an implicit Euler scheme of the system ( |12a[ )-( |T2^ . The discretization hneness is 
given by r := where M e N. We set := := (u^, c°, z°) and construct for 

TO S { 1 ,..., M} recursively by considering the functional 


C, z) := £e{u, c,z)+'R 


m-1' 
^M,e 


T + 


i' 


«:.‘i 


2 

Vo 


2 r' 


m-l||2 

^M.e llL2(n)- 


The set of admissible states for ^ is 


g™ ^ := = (u,c,z) G Wid(L!;K") x H^il) x W^^P{n) 

with u|r = b{mT)\r, / c — c° da; = 0 and 0 < z < z^~^ a.e. in Q 

Jn 

A minimization problem for the functional c, z) = E^ ^{u, c) = d|Vcp+Wch(c)+Wei(e(u), c) da;+ 

^||c — containing a weighted K”))*-scalar product (•, •)l has been considered in |Gar00| . 

However, due to the additional internal variable z, the passage to M —)■ 00 becomes much more involved. 

In the following, we will omit the e-dependence in the notation since e G (0,1] is hxed until Section 

15.31 


Lemma 5.4 The functional E'ffj has a minimizer q'fjj = z"^) G Qffj- 
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Proof. The existence is shown by direct methods in the calculus of variations. We can immediately see 
that is closed with respect to the weak topology in x x W^’^(r2). Furthermore, 

we need to show coercivity and sequentially weakly lower semi-continuity of defined on 


(i) Coercivity. We have the estimate 


Therefore, given a sequence {qk}keN in Qm with the boundedness property £^( 9 ^) < C for all 
fc e N, we obtain the boundedness of Uk in by Poincare’s inequality {uk has fixed boundary 

data on T), the boundedness of Cfc in H^{Q) by Poincare’s inequality (f^Ckdx is conserved) and 
the boundedness of Zk in by also considering the restriction 0 < Zfc < 1 a.e. in 

(ii) Sequentially weakly lower semi-continuity. All terms in E)(} except Wch(c) dec and Wei(e(it), c, 2 :) da: 
are convex and continuous and therefore sequentially weakly l.s.c.. Now let {uk,Ck,Zk) {u,c,z) 
be a weakly converging sequence in In particular, —>■ 2 in LP[Vt), Zk ^ z a.e. in ft and 

Cfc —>■ c in L^{kl) as fc —>■ 00 for all 1 < r < 2* and Ck ^ c a.e. in fl for a subsequence. Lebesgue’s 


generalized convergence theorem yields J^Wchick) dx —>■ J^Wchic) dx using (A 6 I. The remaining 
term can be treated by employing the uniform convexity of Wei(-, c, 2 ) (see (A 1 )): 


/ Wei(e(Mfc), Cfc, 2 fc) - Wei(e(M), c, z) da: 

JQ 

= / Wei(e(it),Cfc,2fc) - Wei(e(M),c, 2)da:-|- / IFei(e(ufc), Cfe, 2fc) - IFei(e(u), Cfc, 2fe) da: 

an an 

> / Wei(e(u),Cfc, 2 fc) - Wei(e(M),c, 2 )da:-|- / 5ekFei(e(M), Cfc, 2fe)(e(ufe) - e(u)) da:. 

an an 

^ "V 

—>•0 by Lebesgue’s gen. conv. theorem and 

The second term also converges to 0 because of deWei{e{u),Ck, Zk) —t deWf,\{e{u), c, 2 ) in L‘^{Q) (by 
Lebesgue’s generalized convergence theorem and (11a)) and e(ufc) — e{u) ^ 0 in L'^{kl). 


Thus there exists q^ = {u^,c^,z^) G Q^j such that E™(g™) = inf^gg™ E™ (g). 


The minimizers q'^ for m £ {0,..., M} are used to construct approximate solutions qM and qM to our 
viscous problem by a piecewise constant and linear interpolation in time, respectively. More precisely. 


qM{t) ■— qMy 

mit) := Pq^ + {1-P)q^-^ 

with t G ((m — l)r, mr] and /3 = retarded function q]^ is set to 


9m (^) 


qM{t-T), if tG[T,T], 
9°, iftG[0,T). 


The functions 6 m and bj^ are analogously defined adapting the notation := birnr). Furthermore, the 
discrete chemical potential is given by (note that dtCM{t) G Vb) 

■■= -(-A)”^ (dtCMit)) + XM{t) (30) 


with the Lagrange multiplier Am originating from mass conservation: 

Xniit) := / dcWch{cM{t)) + dcWei{e{uM{t)), CM{t),ZM{t)) dx. (31) 

Jn 
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The discretization of the time variable t will be expressed by the functions 

dnif) '■= min{mr | m € Nq and mr > t}, 

min{(m — 1)t | to S Nq and tot > t\. 

The following lemma clarifies why the functions <1m 9 m are approximate solutions to our problem. 


Lemma 5.5 (Euler-Lagrange equations and energy estimate) The tuples qm, Qm 9m satisfy 
the following properties: 

(i) for all t G (0, T) and all ^ G 

[ {dtCMit))Cdx = - [ VpMit) -VCdx, (32) 

Jn Jn 

(ii) for all t G (0, T) and all ^ G H^{D,) 

[ h'M{t)Cdx = [ \7cM{t) ■ VC + (9cWch (cm (O)Cdx 
Jfl JQ. 

+ [ dcWei{e{uM{t)),CM{t),ZM{t))C + e{dtCM{t))Cdx, (33) 

JQ 

(Hi) for all t G [0,T] and for all ( G VLpd(f2;K") 

0= f deWei(e(uM(t)),CM(t),ZM(t)) : e(C) + e|V mm(0PVmm(0 : VCdx, (34) 

JQ 

(iv) for all t G (0,T) and all C G IT^’^(n) such that there exists a constant v > 0 with 0 < j^C + ZMit) < 
zf^{t) a.e. in 

0 < / \VzM{t)\P~‘^VzM{t) + d^W^i{e{uM{t)),CM{t),ZM{t))C, - aC, + l3{dtZM{t))C,Ax, (35) 

JQ 

(v) for all t G [0, T] 

rdM{t) pdMit) p ^ ^ 

£6iqM{t))-\- n{dtZM)ds-\- / -|5iCMp + ;^|V^Mpcixds 

Jq Jo Jn ^ ^ 

pdM{t) p 

<4(9°)+/ / deWei{e{uM+ b-bM),c'^,ZM) : e{dtb)dxds 

Jo Jn 

pdM{t) r 

+ e / \\7u]^+\7b-\7b]^\'^\7{u]^+ b-b'lf) :\7dtbdxds. (36) 

Jo JQ 

Proof. Using Lebesgue’s generalized convergence theorem, the mean value theorem of differentiability 
and growth conditions (11a), (A4)-(A6), we obtain the variational derivatives of 4 with respect to u, c 
and z: 


(44(g), C) = 

f deWei{e{u), c, z) : e(C) +e : VCdx for C G 

JQ 

(37a) 

( 44 ( 9 ), C) = 

[ Vc-VC + 44ch(c)C + (9cWei(e(M),c,z)Cda; for C G i?^(U), 

Jq 

(37b) 

(44(g), C) = 

[ \\7z\P-^\/z ■ VC + d,Wei{e{u),c, z)C dx for C G W^’P{n). 

JQ 

(37c) 


To (i)-(v): 
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(i) This follows from (30). 

(ii) fulfills (c1cEm(9m): Ci) = 0 for all S Vb and all m S { 1 ,..., M}. Therefore, 

0 = (dcfe(qM(i)),Cl) + {dtCM{t),Cl)v^ +£( 9 tCM(<).Cl)L 2 (n)- 


On the one hand, definition (30) implies 


{dtCM{t), Cl)Vo = ((“^) ^ i.dtCM{t))Xi)L^(n) 

= { — + XM(t), Cl)L2(n) 

= Ci)l2 (0) 

and consequently 

0 = {dc£e{qM{t)), Cl) - M{t), Cl)+ e{dtCM{t), Cl) L^{n) for all Ci G Vb- 

On the other hand, definitions (30) and © yield for C 2 = C with constant (7 G M: 

(dc£e(gM(f)), C2) — (MM(f), C2)i2(0) + e{dtCM{t),C2)L^{Q) 

= C'£"( 0 )AM(f) + ((—A) ^ {dtCMit )), C2)L'^ip.) — (AmW, C2)L2(n) + 0 


(38) 


=0 




= 0 . 


(39) 


Setting = C — J^C and C2 = / C) inserting (37b) into (38) and (39), and adding (38) to (39) shows 
finally (ii) (cf. [GarOOl Lemma 3.2]). 


(hi) This property follows from (37a) and 0 = (d„E)(J^(( 7 ](}), ^) = (du£e(qM),C) for all C G Wj!;’'^(0; K"). 
(iv) By construction, zJJ) minimizes E'^{u'^, d^, •) in the space W^’P{fl) with the constraints 0 < 2 ; and 


m—1 


< 0 a.e. in O. This implies 


— {^z£siqM),C ~ ^m) ~ ( 


^m — 1 


X-zZ) <0 

L 2 (n) 


(40) 


for all C G IT^’^(fl) with 0 < C < z’^ ^ a.e. in O. Now, let the functions C G and v > 0 


with 0 < + ZM{t) < Zj^^{t) a.e. in Q be given. Since 1 / > 0, we obtain from (40): 

-{dz£eiqM(t))X{t)) - {din{dtZM{t)) ,C{t))L^n) < 0 . 


This and (37c) gives (iv). 


(v) Testing EJJJ with q = ^ 


M ’ '-M ’ 


^) and using the chain rule yields: 


m-l 


r + y^Wcm - c^-%^ + - XXTlHP 


4(9M)+7^ 

\. / y 

— f -U f _L _ F 

— ^sKHM ) ^ ^€\^M ^ ’'^M / ) 

d 


m — 1 I j-,f ^Tn—1 


bis) - h' 


M ’ ''M ’ 


) ds 


prri' 

= 4(gM-^)+/ 

J{m-1)T OS 

= 4(9m-') 

piriT n 

+ / / deWeMu^-^ + bis) - b^-^),clf\zZ-^) : eidtb) dxds 

J (m—l)r J Q 
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prriT p 

+ e / + V&(s) - + b{s) - b^-^) : V5t6dxds. 

*/ (m— l)r J Q 


' (m— l)r ^ Q 

Summing this inequality for k = 1,... ,m one gets: 


£e{qZ) + Y.r[n 


— z^~^ \ 1 


k=l 




Vo 


rk _ 


L2(n)y 


pmr p 

<^e( 9 °)+/ / 9eVhei(e('u^+ 6 -: e(9t&)da;ds 

Jo Jf 2 

^mr ^ 

+ £ / / |Vu^^ + V 6 - V5^PV(m^ + 6 - 6 ^) : V9(6dxds. 

Jo Jn 


Because of 


/o Jn 

,fc ||2 




= II^/^mIIl 2 (q) by (30), above estimate shows (v). 


The discrete energy inequality ( |36| ) gives rise to a-priori estimates for the approximate solutions. 

Lemma 5.6 (Energy boundedness) There exists a constant C > 0 independent of M, t and £ such 
that 


rdM{t) rdM^t) p i 

4(<?M(i))+/ n{dtZM)ds+ / -|dtCMp + xlV^Mpdsds < C(£:e(g°)-f 1). 

Proof. Exploiting (A3) yields the estimate (C > 0 denotes a context-dependent constant independent 
of M, t and e): 

/ deWei(e(u^(s) -f b{s) - blj{s)), cfjis), z^{s)) : e{dtb{s)) da; 

Jn 

< C'||V5t6(s)||ioo(n) f Wei{e{ulj{s)),clj{s),z^{s)) + \e{b{s)-blj{s))\ + ldx. (41) 
Jn 


rdM{t) 


In addition, 


[ |Vmm('S) +^Ks) - V 6 jv,^(s)pV(uj\^(s) -f b{s) - bjyf{s)) : \'dtb{s)dx 
Jn 

<q|Vd* 6 (s)|U»(o) [ \Vui,is)\^ + \V{bis)-b-Mis))\^dx. 


(42) 


To simplify the notation, we define the function: 


._ \ £e{qM{t)) + TZ{dtZM)ds + §|dtCMp + 5 |V/rMpda;ds, ifte[ 0 ,r] 


pdM(t) 


7(0 := 


[£e{q°), 


ate [—r, 0). 


Using ( |41[ ) and ( |42[ ), the discrete energy inequality p 6 | ) can be estimated as follows: 

ndM{t) 


paM{t) 

lit) < £eiq°) + C / ||Vat6(s)||icx,(Q)4(g- (s))ds 

Jo 

-I- C'llII|V(5 — bj^,j)\ + |e(&— bj^)\ + l||ioo([Q 
< £M) + CJ II va*6(s + r)|U.o(n)4('?M(s)) ds + C 
<£M) + C f ||Vdt&(s + r)|U=o(n)7(s)ds + C. 
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Gronwall’s inequality shows for all t G [0,T] 


7 (t)<C + £,(q°) + C J' (C + £:e(g°))||V5t6(s + T)|U=e(a)exp(^^‘||V9t6a + r)|U=o(o)dZ^ ds 

<c(4(a + i). 


Corollary 5.7 (A-priori estimates) 

fi) l|uM||L“([0,T];VKi.^(n;R’*)) < C, 
(a) ||cM||L°°([O.T];ffi(n)) < C, 

(Hi) ||2M||L“([0,T];H^i.P(n)) < C, 


There exists a constant C > 0 independent of M such that 
(iv) \\dtCM\\L2{nT) ^ 

(v) ||'9t^M||L2(Q^) < C, 

(vi) ||MM||L2([o,r];ffi(a)) < C 


all M gN. 


for 


Proof. We use Lemma 5.6 


The boundedness of {V{uMit) — bAiit))} in L^{n;] 


) and UMit) — bM(t) G 
by Poincare’s inequality. The boundedness of {VcM(i)} and mass conser- 


yield (i) 

vation imply (ii) by Poincare’s inequality. The boundedness of {'VzM{t)} in L^’(O) and 0 < ZM{t) < 1 
a.e. in ft for all M and all t G [0,T] show (iii). The properties (iv) and (v) follow immediately from 
The boundedness of {V/tm} in and qiM(t) dx} with respect to M and t show 


Lemma 


5.6 


(vi) by Poincare’s inequality. Indeed, HMif) dx} is bounded with respect to M and t because of (33) 
and (32) tested with C = 1- 


Due to the a-priori estimates we can select weakly (weakly-*) convergent subsequences (see Lemma 


5.8). Furthermore, exploiting the Euler-Lagrange equations of the approximate solutions, we even attain 


strong convergence properties (see Lemma 5.9 and Lemma 5.11). 


Lemma 5.8 (Weak convergence of the approximate solutions) There exists a subse¬ 
quence {Mk} and elements {u,c,z) = q G and fi G L^([0, T]; iJ^(n)) with c(0) = c°, z(0) = z^, 
0 < z < 1 and dtz <0 a.e. in fix such that the following properties are satisfied: 


(i) ZM,,z-^^^ z in L°°{[Q,T]-W^^^{n)), 
ZMt,{t),z'fj^{f) z{t) in IT^’^’(Sl) a.e. t, 

ZM,^z in H\[Q,T]-,L\9.)), 

(ii) UM, ^uin L°°([0,T];lFi’4(n)), 
as k ^ oo. 


(iii) m L°°([0,r];iLi(fl)), 

c{t) in H^{n) a.e. t, 
“t c a.e. in fix and 
CMk c in H^{[0,T];L‘^{n)), 

(iv) h'Mk ^ H in L‘^{[0,T]-H^{n)) 


Proof. To simplify notation we omit the index k in the proof. 

(iii) Since cm is bounded in L^{[0,T]; H^{fl)) and dtCM is bounded in L^(Jlx) we obtain cm —t c in 
Lf{Vlx) as M —7> oo for a subsequence by a compactness result from J. P. Aubin and J. L. Lions 
(see |Sim86] L Therefore, we can extract a subsequence such that CM{t) —t c{t) in L^(n) for a.e. 
t G [0,T] and cm ^ c a.e. on ^x- We denote this subsequence also with {cm}- The boundedness 
of {cM(i)}MeN in H^{Tf) even shows CM{t) c{t) in for a.e. t G [0,T]. In addition, the 

boundedness of {cm} in -b°°([0, T]; iL^(f2)) shows cm ^ c in L°°([0, T]; iL^(n)). Furthermore, we 
obtain from the boundedness of {SjCm} in L^(11 t) for every t G [0,r]: 

||cM(i) - CM(l)||Li(a) = \\cM{dMit)) - CM(i)||Li(n) 

dM{t) 

ll' 9 tCM(s)||Li(n) ds 

< C{dM{t) - t)'^^'^\\dtCM\\mQT) ^ 0 as M oo. 
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Lebesgue’s convergence theorem yields \\cm — CM\\L^{aT) —>■ 0 as M —)• oo. Analogously, we obtain 
IICm — c^j||z,i(nT) —S' 0 as M —>■ oo. Thus, the convergence properties for cm also holds for cm and 
c^ with the same limit c = c“ = c a.e. . The boundedness of {cm} in L‘^{n)) shows 

Cm ^ c in 7J^([0,T];for a subsequence. 

(i) We obtain the convergence properties for {zm} with the same argumentation as in (iii). Note that 
the limit function is also monotonically decreasing with respect to t. 

(ii) This property follows from the boundedness of {um} in L°°([0,T];K")). 

(iv) This property follows from the boundedness of {hm} in L‘^{[0,T]; H^{n)). ■ 


In the sequel, we take advantage from the elementary inequality (x, y are elements of an inner product 
space X with scalar product (•,•)) 

Cnc\\x-y\\‘^ < ((||x||«"^x- \\y\\'^~^y),x-y) (43) 

for a constant Cue > 0 depending on X and q>2. To see this, ( |4^ is equivalent to 
Cue < {b, ||a + 6||'^“^(a + 6) — ||a||'^“^a) for all a,b G X, ||6|| = 1 
by introducing the variables a := x/||x — y\\ and b := (x — y)/\\x — y\\ for x ^ y. This is equivalent to 

Cue < I|a + 6||'^“^ + (&, o) (||a + b\\‘^~'^ - ||a||'’“^) for all a,bG X, ||6|| = 1. (44) 

Now the equivalence ||a + 5|| < ||a|| (a,5) < —^H^IP gives the estimate: 


|a + b\r^ + (6,a) (||a + b\r^ - Haf-^) > ||a + b\r^ + -||&f (Hair-^ - ||a + b\r^) 

= l\\a+br-^ + l\\ar-^ 


Since ||&|| = 1, the right hand side is bounded from below by a positive constant and therefore (44) 
follows. 


Lemma 5.9 There exists a subsequence {Mk} such that xl in L'^([0, T]; M")) as 

k —)■ 00 . 


Proof. We omit the index k in the proof. 

Applying ( |Al[ ) , taking inequality (431 for <7 = 4 into account and considering (341 with the test-function 
C = UM{t) - u{t) - bM{t) -f b{t), we get 

V\\e{UM) — e(u)||i 2 (Q,j,.]I{r.Xri) -f eCuc||VUM — VM||}^4(Q,j,.]I{r.Xri) 

< / (l9eWel(e(MM),CM,ZM) - 9elTel(e(lt),CM,-ZM)) : (e(uM) - e(u))dxdf 

J 

+ £ (iVuMpVuM — iVupVu) : {Xum — Vm) dxdt 

J Ot-- 

= / 5eWei(e(uM),CM,ZM) : e(C)-f elVuMpVuM : VCdxdt 


=0 by l|34|l 


-I- / deWe\{e{uM),CM,ZM) ■ {e{bM) - e{b))dxdt+e | VwmPVmm : (V 6 m - V 6 ) dxdt 

J Q J- J Q-j” (**) 


(*) 


(**) 
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/ {df.Wf.\{e{u),CM,ZM) ■ {e.{uM) — e{u))AxAt—e / | : (Vum — Vu) dxdi . (45) 

J CIt J r2x' 






Since deWe\{e{uM), cm, zm) is bounded in L^(riT;R"^”) (by (11a) and Corollary |5.7| as well as e(&M) —> 
e(b) in we obtain (*) —>■ 0 as M —)■ oo. The boundedness of |VumPVum in M"^") 

by Corollary 5.7 and V 6 m —t V& in L‘^(11 t; M"^") lead to (-a^) —>■ 0. We also have deWe\{e{u), cm, zm) —t 
deWei{e{u), c, z) in L^(nj’; K”^”) by ( lla| ) and Lebes gue ’s generalized convergence theorem. Further¬ 
more, e{uM) e(u) in L‘^ {VIt ', R” x K”) by Lemma 5.8 This gives (* * *) —>• 0. Since Vum ^ Vw 
in L"‘(flr;K") by Lemma 


5.8 


we obtain ( 7 k-***) —>• 0. Therefore, (45) implies e{uM) —t e(it) in 
) and Vwm —t Vu in L"‘(Ot; as M 00 . Poincare’s inequality finally shows um —t u 

in L"‘([0,T]; lF^’'’‘(fl;M”)). Now, we choose a subsequence such that UM{t) —t u{t) in 1F^’"‘(0; M") for a.e. 
t G [0,T] and um —t u a.e. in O 7 ’. We also denote this subsequence with {um}- 

Analogously, we obtain a u~ G L'^{[0,T]-,W^’‘^{il.)) satisfying -G u~ with the same convergence 
properties. We will show u = u~ a.e. . Consider (34) for qMit) and for 


0 = 


df,We\{e{uM),CM,ZM) : e(C) + e|V mmPVum : VCdxdt, 

J Qj’ 

0 = / deWe\{e{u~^),c1i,z~^)-e{C,)+£\Vuli\^VuM-.VCdxAt. 
J Qt- 


(46a) 

(46b) 


We choose the test-function C,{t) = UM(t) — Uf^{t) — bM{t) + bf^{t) G VFp’‘^(f2). An estimate similar to 
(45) gives: 

v\\e{uM) — e(M^)||i 2 (Q,j,) -|- eC-^^^\\VuM — 

< / {deWe\{e{uM), CM, Zm) “ 3elTel(e(lt)(p), Cm, Zm)) ■ {e{uM) - c(u^)) dxdt 

J 

+ e (iVuMpVuM - I'^ulfl'^Vu]^) : {Vum - Vu)(p) dxdt 

J Ot-- 

= / 5elTei(e(MM),CM,ZM) : e(C)+ e|VMMpVuM : VCdxdt 


=0 by (|46a| 


/ deWe\{e{uj^^),Cfj,Zj^) : e(C) + £| : VC dxdt 

J Q J- 


=0 by ( |46b[ ) 

-I- / (5eWei(e(M)(p),c)(p,z)(p) - 9eWei(e(u)(p), CM, ^m)) : (e(uM) - e(u)(p)) dxdt 

J 0^1 

-I- / (5eWei(e(MM),CM,2:M) “ deWei{e{u~^), c]^, z~^)) : {e{bM) - e{blf))dxdt 

J 

+ e f (iVuMpVuM - |Vu)(ppVu)(p) : (V6 m - V6 )(p) dxdt. 

J 

Observe that deWei{e{u~^), c~j^, z'^) — deWei{e{u~^),CM, zm) —>■ 0 in L‘^{Qt) by Lebesgue’s generalized 
convergence theorem (using growth condition (11a), Lemma 5.8 and convergence properties of um and 
u]^) as well as e{bM) — c{b~^) —)■ 0 in L'^iTlT', K"^”) and V6 m — 0 in L‘^{nT', M"^”). Hence, each 

term on the right hand side converges to 0 as M —)■ 00 I 


Lemma 5.10 There exists a subsequence {M^.} such that Cm^,Cj^,i^ -G c in L^([0, T]; 77^(0)) as k ^ 00 . 
Proof. We omit the index k in the proof. 

implies CM{t) -G c(t) in L^*/^+^(0) for a.e. t G [0,T]. Using Corollary 5.7 and Lebesgue’s 


Lemma 


5.8 
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convergence theorem, we get Cm —t c in Next, we test (33) with ^ = CM{t) and integrate 

from t = 0 to t = T. Then we use Lebesgue’s generalized convergence theorem, growth conditions (A4) 
and (A6) as well as Lemma 5.8 to obtain 

/ |VcMp dsdt —>■ — / dcWc\i{c)c + dcWe\{e{u),c,z)c + e{dtc)c — iJLcAxdt 

J O'j' J Ot-- 

as M —>■ oo. On the other hand, we test (33) with c{t) and integrate from t = Q to t = T. Note 
that c S {VLt) and dcWc\i{cM) —t dc^chic) in LF ‘*^as M —)• oo by Lebesgue’s generalized 
convergence theorem. Hence, we derive for M —> oo: 

/ |Vcpda:dt = — / dcWc\i{c)c +dcWf.\{e{u),c,z)c +£{dtc)c —^cdxdt. 

J J 

Therefore, cm -t cin L‘^{[0,T]-H^{n)) as M oo. The convergence \\cm\\l^{[o,t]-,h^(Q)) -t ||c||L 2 ([o.T];i?i(f 2 )) 
implies \\c^\\Lmo,T]-min)) ||c||L 2 ([o,T];//i(n))- We also have ^ c in L'^{[Q,T]-, H^{n)) (by Lemma 

(ii)) and consequently c'^j —)■ c in L^([0, T]; iJ^(r2)) as M —>• oo. 


5.8 


Note that in connection with Corollary 5.7 we even get for each g > 1 

CM,clj^cm L'^i[0,T];H\n)) 

for a subsequence as M ^ oo. 

Lemma 5.11 There exists a subsequence {Mk} sueh that ZMk^ z in U‘{[t),T]\W^’'P{Vl)) ask^oo. 

Proof. To simpli fy no tation we omit the index k in the proof. 

Applying Lemma 5.2 with / = C = and /m = gives a sequence of approximations {Cm}mgn C 
L^’([0 , T]; 1T_J:’^(H)) n L°°{Qx) with the properties (note that we have z^{t) z{t) in W^''p{VL) for a.e. 
t € [0,T] by Lemma: 


(m^z in LP([0,T];W^’P(n)) as M ^ oo 
0 < Cm < -^M M e N. 


(47a) 

(47b) 


We test (35) with C = Cm(0 — ZM(t) for i/ = 1 (possible due to (47b)), integrate from t = 0 to t = T and 
use (43) to obtain the following estimate: 


Cue / IVzm — Vz|^ dxdt 

[ (IVzmT’^Vzm - |Vz|P-2Vz) • V(2 m - z) dxdt 

[ IVzmT'^Vzm • V(zm - CM)da;dt 

J Q'j-' 

+ f IVzmT’^Vzm • V(Cm -z)-\Vz\p-^Vz ■ V(zm - z) dxdt 
J 

; / (9^H4i(e(uM),CM,2:M) - Q! + ,59tZM)(CM - ■ZM)da:dt 

J Q'j’ 

+ f IVzmT’^Vzm ■ V(Cm -z)-\Vz\p-'^Vz ■ V(zm - z) dxdt 
< \\dzWei{e{uM),CM, zm) — a + f3dtZM\\L^{nT) IICm — ZM\\L^{nT) 

bounded by fAst and Cor. EH 

+ llVzM||M(n,j,) I|VCm - Vz||M(aT) - / • V(zm - d dxdt. 

^ ^ ^ J VLt 

bounded by Cor. |5.7| 


< 

< 


< 
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Observe that V^m — Vz —>■ 0 in L 2 ’( 07 ’;K") and Cm — zm —>■ 0 in (by property (47a) and by 

Lemma 5.8) as well as Vzm ~ Vz ^ 0 in L^(Or;IR") by Lemma 5.8 Using these properties, each term 
on the right hand side converges to 0 as M —)• oo. 

We also obtain 

\\zm\\lp{[o,t]-w^’P{Q)) —t lkllLp([o.T];rvi.P(n)) fr om | |zM||Lp(fo,ri:tvcp(o)) —t ||z:||LP([o,T];Wi.p(n)). Because 
of z^ ^ z in LP{[0,T];W^’P{n)) (by Lemma 5.8 (i)) we even have z^ —)■ z in LP{[0,T];W^’^{n)) as 
M —)■ 00 . I 


In conclusion, Corollary |5.7[ Lemma 5.8 Lemma 5.9 Lemma 5.10 and Lemma 5.11 imply the following 
convergence properties: 

Corollary 5.12 There exists subsequenee {M^} and an element (u,c,z) = q G Q'' with c(0) = c° and 
z(0) = z° such that 

(Hi) -Gu m L‘^([0,T]; 1U^’'^(0;K”)), 

UMk{t),u^^{t) -G u{t) in lUi’'‘(0;]R"') 
a.e. t, 

UMk , ^Affc 


(i) ^ z in LP{[0,T]-,W^’Pm, 

^M,(t),2^,(0 ^ z{t) in W^’P{n) a.e. 


^Mk 


^Mk 

{t),z 

i, 


^Mk 


^Mk 



(in) ^lMk ^ H m L^([0,r];i7i(O)) 

(v) dcWc\,{cMk) -t 9cWch(c) in L'^{VtT) 


z a.e. in VLt and 

(a) CMk,CM^ -G c in L2*([0,T];i7i(O)), 

CMk{t),c~j^^{t) -G c(t) in H^{Vl) a.e. t, 

CMk^c]^^ —>■ c a.e. in fix and 
CMk ^ c in H^{[0,T];L‘^{^)) 
as k ^ oo. 

The above convergence properties allow us to establish an energy estimate, which is in an asymptotic 
sense stronger than the one in Lemma |5.5| (v). We emphasize that (36) has in comparison with (481 no 
factor 1/2 in front of the terms P\dtZMr7s\dtCM\'^ and |V/rMp. 

Lemma 5.13 (Precise energy inequality) For every 0 < ti < t 2 < T: 

ndM (* 2 ) 


paM\^2) p 

Se{qM{t2))+ / / -oStZM +/?|9tZMp+e|9tCMp + iV^Mpdxds - 5E(g)/f(tl)) 

Jn 

pdM{t2) p 

< / i9eWei(e(M/;f+ 6-: e(at6 )dxds 

dd^iti) Jn 

pdM{t2) p 

+ e / \\7ulf + \7b-\7blj\‘^\7{u~^ + b-blj):\7dtbdxds + KM 

''d.fjdi) JQ. 


(48) 


with Km ^ d as M ^ 00 . 

Proof. We know W^{q'^) < ~ z]^). The regularity properties of the functions 

6, Cm and zm ensure that the chain rule can be applied and the following integral terms are well defined: 

£M,cZ,z^) 

<£,{uZ-^ + bZ-bZ-\cZ,zZ) 

= £ZuZ-\cZ-\zZ-^) 

f M h'^ — _ f 

I r — ^ I 1\ c I 7m L’Ti—1 ^m—1 ^m—l\ 

+ ££{nM d-bM — bj^ ,Cm, Zm) — £e[UM +bM~bM ^Cm,Zm ) 

pmr 

+ / {dMuZ-^+b{s)-bZ-\cZ-\zZ-^),dtb{s))^H^y^H^ds 

J (m— 1 )t 
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+ / {dMuZ-^+b^-bZ-\cM{s),zZ-^),dtCMis))^my^mds 

J (m— 1 )t 

pmT 

+ / {dz^e{u^ ^ ~ C^: ^m(s)), 9tZM(s))(vi/i,p)*xlVi.P ds. 

J (m— 1 )t 


Summing from m = + 1 to yields: 

4(<7m(<2)) - SeigMih)) 

pdMit2) r 

<e / |V(u^ + &-5^)|^V(w^ + 5-6^) : V9t5da;ds 

rdM{t2) r 

+ / 5eWei(e(w^ + 6-&^^),c^,0^) : e(at6)dxds 

pdM{t2) p 

+ / 5cWei(e(uM + &M - ^M)>dM,z^)9(CMda;ds 

Jn 


rdM (^ 2 ) 


j(ti) Jq 


(*) 


Vcm ■ '^dtCM + dcWch{cM)dtCM dxds 


(★★) 


pdM{t2) p 

+ / 5zWei(e(M^+ 6 m - &M)>CM,-SM)9tZM + IVzmT'^Vzm • VStZMdxds. (49) 

''d-Mdl) 

" -V-' 

(*★*) 

By using convexity of a; 1 —> we obtain the following elementary inequality 

i\VzMit,x)\P~'^VzM{t,x) - \VzMit,x)\P~'^VzMit,x)) ■ VdtZM{t,x) < 0 . 


This estimate and (35), tested with C = —dtZM^t) for v = t and integrated from t = 0 to t = T, lead to 
the estimate: 


pdM{i2) r 

(***)< - / -adtZM + PldtZMl"^ dxds 

dd-{ti) Jn 


Jn 

f'dM{i2) 


paM{^ 2 ) p 

/ / idzWel{e{uJ,i +bM - bJf ),CM,ZM) - dz;Wel{e{uM),CM,ZM))dtZMdxds . 

Jd^di) J^ 


Furthermore, 


pdM{t 2 ) p 

W< /_ / - 

•'dM^i) dn 


5cWei(e(uM), CM, ZM)dtCM dxds 


pdM{t 2 ) p 

+ / {dcWel{e{uJf + bM - blj),CM,ZM) - 9cWe\{e{uM),CM,ZM))dtCMdxds . 

Jdj^^^{ti) Jn 
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Using the elementary estimate {Vcm — ^cm)^ dtCM < 0, we obtain 

pdM{t2) p 




^(ii) JQ. 
pdM{t2) 


Vcm ■ ydtCM + dcWchicM)dtCM dxds 


'(ijvrbi) 


{dcWch{cM) - dcWch{cM))dtCM dxds . 


Hence, applying equations (33) and (32) shows 

pdM (* 2 ) 




{dcSe{qM), dtCM){H^)*xH^ ds = 


pdM{t2) 

'dj^{ti) do, 
pdMit2) 


- £|9tCMp dxds 


paMV'2) r 

/ / -|V^Mp-el^tCMpdxds. 


Thus, 


pdM{t2) r 

(*) + (**)< / / -|V/rMp - £|9tCMpdxds ++ Km- 

Jfl 


Lebesgue’s generalized convergence theorem, the growth conditions ( |A4[ ), (A5), (A6) and Corollary 5.12 
ensure that k\,j, k\^ and k|^ converge to 0 as M —00. Here, we want to emphasize that we need 
boundedness of dtCM and dtZM in LF‘{0.t) and the convergence e(uM) —t e(u) in which we have 

only due to the regularization for every fixed e > 0 as M —>■ 00 (see Corollary 5.12). To finish the proof, 
set Km ■= + k\j. ■ 

We are now in the position to prove the existence theorem for the viscous case. 

Proof of Theorem |4.4[ The proof is divided into several steps: 


(i) Using growth conditions ( |A4[ ), (A6), (Hal, Corollary 5.12 and Lebesgue’s generalized convergence 
theorem, we can pass to M —> 00 in the time integrated version of the integral equations (32), (33) 
and (34). This shows (i) and (ii) of Definition |4.3| 


(ii) Let 0 < < ^2 < r be arbitrary. Because of d^(ti) < U < H < dM(^ 2 ), Lemma 5.13 particularly 

implies 


£s{qMit2)) 


-adtZM + PldtZM^ +£|i9tCM|^ + |V^m|^ dxdt - 


rdM(t2) p 

/ / 9eWei(e(u((f + 5 - c~^,zm) : e{dtb) dxdt 

dd-^{ti) Jn 

pdM{t 2 ) p 

+ e / + \7b-Vbljf\7{ulj+ b-b^j) :\7dtbdxdt +Km 

dd^(ti) Jo, 


< 


(50) 


with Km —>■ 0 as M —>• 00 . Due to growth condition (A2|, (A6), Corollary 5.12 and Lebesgue’s 
generalized convergence theorem we obtain 


£e{qM{t)) £eiq{t)) and £e{qM{t)) £e{q{t)) 


(51) 


as M —)■ 00 for a.e. f S [0,T]. A sequentially weakly lower semi-continuity argument based on 
Corollary |5.12| shows: 

liminf [ [ -adtZM + P\dtZM\'^ + e\dtCM\'^ + \'^hm\^ dxdt 

M^qo 7j-2 
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> / a{z{ti) - z{t 2 )) Ax + f f /3|i9tzp + e|9(Cp + |V/xp da;di. (52) 

Jo, J ti JQ. 

Growth condition ( |lla ), Corollary 5.12 and Lebesgue’s generalized convergence theorem show: 

deWe\{e{ulj + 6 - zm) ^ deWe\{e{u), c, z) in L°°{[0, T]-,L'^{n)), 

\Vum + Vb- V6)(^pV(w)(f + 6 - &m) ^ \Vu\^Vu in L~([0, T]; 

Since e(dtb) G T]; L^(r2)) and Vdtb G T]; L'*(n)) we get: 

pdM{t 2 ) r 

/ / deWe\{,e{u'^i + b- blj), c~^, zm) : e{dtb) dxdt 

->• f I deWei{e{u),c,z) : e{dtb) dxdt, 

Jti Jq 

pdM{t2) p 

/ / \^ulj+Wb-Wb~M\'^V{ulj + b-blj)-.Vdtbdxdt 

Jq 


(■*2 


-)> 


f JQ, 


|Vm|^Vu : Wdtbdxdt. 


(53) 


Now, using (51), (52) and (53) gives (iv) of Definition 4.3 by passing to M —>• oo in (50) for a 
subsequence. 

(hi) Let ( G LP([0, r]; VLI’^(D)) n be a test-function with {^ = 0} 3 {z = 0}. Applying 

with f = z and Jm = zm and C, = —C gives a sequence of approximations {Cm}m6N C 


Lemma 


5.2 


LP{[0,T]-,W^’^{n)) n L°°(Dr) with the properties: 

Cm ^-C in LP{[Q,T];W^'P{n)) as M ^ oo, 

0 < i'M,tCM{t) < ZM{t) a.e. in D for a.e. t G [0,r] and all M G N. 


(54a) 

(54b) 


Let Cm denote the function —Cm- Then, (54b) in particular implies 0 < i'M,tCM{t) + ZM{t) < Zj^{t) 
a.e. in D for a.e. t G [0,T]. Now, (35) holds for C = Cuit)- Integration from t = 0 to t = T 
and using growth condition (A5), Corollary 5.12 and Lebesgue’s generalized convergence theorem 
as well as the strong convergence (54a) yield for M —>■ oo: 


f |Vz|^ • VC +i92lTei(e('u),c, z)C —Q!C + /3(i?t-z)Cda:dt < 0. 

J 


(55) 


(iv) Property (55) implies that 

- [ |Vz(<)|P"^Vz(t) • VC -b {d^Wei{e{u{t)), c{t), z(t)) -a + l3{dtz{t)))C,dx < 0 
Jn 

holds for all C G W]_'P{tX) with {C = 0} A {z(f) = 0} and for a.e. t G [0,T]. Applying Lemma 
with / = |Vz(<)|P“^Vz(t) and g = dzWe\{e{u{t)),c{t), z{t)) — a + j3{dtz{t)) shows 

[ |Vz(t)|^“^Vz(<) • VC-b {d:,We\{e{u{t)),c{t),z{t)) - a + /3{dtz{t)))Cdx 


5.3 


> 


> 


'{2(t)=0} 


[d,,Weiie{u{t)), c{t),z{t)) -a + P{dtz{t))]+C dx 


/{2(t)=0} 


[d,,Wei{e{u{t)), c{t),z{t))]+C dx 


(56) 
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for all C € W]^^{Vl). Setting 


'T ■■= -X{z=o}[dzWei{e{u),c,z)\~ 


we get (24) from (561 by integration from t = 0 to t = T and we also have 

(r(t), C - z{t)) = - [ [d^Wei{e{u{t)),c{t), z(t))] + (C - z{t)) dx < 0 

J{zit)=0} 

for any ( £ and a.e. t £ [0,T]. Therefore, (25) is shown. 


5.3 Vanishing viscosity: e \ 0 


For each £ G (0,1], we denote with = {ue,Ce,Ze) £ Q'' a viscous solution according to Theorem 4.4 


Whenever we refer to the equations and inequalities (211-( 27) of Definition 4.3 the variables q = (u, c, z), 
fi and r should be replaced by q^ = {ug,Cg, Zg)^ /i^ and r^. By the use of Lemma 5.14 Lemma |5.15| and 
Lemma |5.16| below, we identify a suitable subsequence where we can pass to the limit. 


Lemma 5.14 (A-priori estimates) There exists a C > 0 independent of e > 0 such that 
(i) ||Ue||ioo([o,T];Hi(n;R")) < C, 


(a) £^/‘*||Me||L=o([0,T];Wi.4(n;R")) < C, 
(Hi) ||c£||i;,oo([o^T];ffi(n)) < C, 

(iy) lke||L“([ 0 ,T];Wi.p(a)) < C, 
for all e £ (0,1]. 


(v) ||9tCe||L2([Q_7.].(^i(Q)).) < C, 

(vi) e^/‘^\\dtCe\\L^nT) < C, 

(vii) \\dtZe\\L^Q.T) < C, 

(via) ||Me||L2([o,T];ffi(n)) < C 


Proof. According to Lemma 5.6 the discretization qM,e of qe fulfills 

pdM (t) 


fdMit) 


£s{qM,s(t)) 


'R-{dtZM,e) ds 


||9tCM,£p + ^|V/XM,6pda;ds < C{Ee{ql) + 1), 


(57) 


where C is independent of M,t,e. By the minimizing property of we also obtain Ss{q).) < £e{qi) ^ 

£i{qi) for all £ £ (0, Ij. Therefore, the left hand side of (57) is bounded with respect to M £ N, t £ [0,T] 
and £ G (0,1]. This leads to the boundedness of 

£e{qs{t))+ f TZ{dtZe)ds+ f / ||i9tCEp + ijV/iepdxds < C7 (58) 

for a.e. t £ [0,T] and for all e £ (0,1]. We immediately obtain (iv), (vi) and (vii). Due to J c^{t)dx = 
const and the boundedness of 1| Vc£(t)llz,^), Poincare’s inequality yields (iii). In addition, (ii) follows 
from Poincare’s inequality. Now, using ( |58[), g rowth conditions (11b) and Korn’s inequality, we attain 
the desired a-priori estimate (i). Due to |22[) and (21) we obtain boundedness of p,s{t) dx. Since 
||V/re(t)|h 2 (Q,j,) is also bounded, Poincare’s inequality yields (viii). 

Finally, we know from the boundedness of {V^e} in L^(Dt) that {dtCg} is also bounded in L^([0, T]; (i7^(D))*) 
with respect to e by using equation (21). Therefore, (v) holds. ■ 


Lemma 5.15 (Weak convergence of viscous solutions) There exists a subsequence {sk} (which is 
also denoted by e) and elements (u,c,z) = q £ Q and p, £ L^([0, T]; 77^(D)) with z(0) = z^, 0 < z < 1 
and dtz < 0 a.e. in £It such that 
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(i) Ze ^ z in L°°{[0,T\,W^’P{n)), 
Ze{t) z{t) in W^’P{Vl) a.e. t, 

Ze ^ z a.e. in Ht and 
Ze ^ z in iJi([0,r];L2(ri)), 

(ii) Me ^ u m L°°([0,r];i7i(O;K")), 
as e \ 0. 

Proof. 


(in) Ce ^ c in L°^{[0,T\,H^{Vt)), 

Ce(t) c(t) in H^{n) a.e. t and 
Ce ^ c a.e. in fir, 

(iv) /ie ^ ^ m L^([0,T];iJ^(fi)) 


(i) This property follows from the boundedness of {z^} in L°°([0, T]; IT^’^(Sl)) and in 

T]; L^(n)) (see proof of Lemma 5.14). The function z obtained in this way is monotoni- 
cally decreasing with respect to t, i.e. dtz < 0 a.e. in fiy. 

(ii) This property follows from the boundedness of {me} in L°°([0,T];iL^(r2;K"’)). 


(hi) Properties (iii) and (v) of Lemma 5.14 show that Cg converges strongly to an element c in 

as e \ 0 for a subsequence by a compactness result due to J. P. Aubin and J. L. Lions (see 
|Sim86| ). This allows us to extract a further subsequence such that Ce(f) —?► c{t) in L^(r2) for a.e. 
t G [0,r]. Taking also the boundedness of {cg} in L°°{[0,T]-,H^{n)) into account, we obtain a 
subsequence with Ce{t) c(t) in H^{n) for a.e. t G [0,r] and — >■ c a.e. in fix as well as ^ c 
in L-([0,r];iLi(fl)). 

(iv) This property follows from the boundedness of {/Tg} in L'^{[0,T]\H^{Q)). ■ 


Lemma 5.16 (Strong convergence of viscous solutions) The following convergence properties are 
satisfied for a subsequence e \0: 


(i) Mg ^ M m L2([o,T];iji(fl;K")), (in) ^ z m LP([0, T]; Ipi’P(fl)). 

(ii) Ce ^ c in L^{[0,T]-,H^(fil)), 

Proof. 


(i) We consider an approximation sequence {M 5 } 5 g(o,i] C L^([0, T]; IP^’^(n)) with 

us ^ u in L^{[0,T]-, H\Q)) as <5 \ 0, (59a) 

us-bG for all (5 > 0. (59b) 


Since e and <5 are independent, we consider a sequence {^£}gG(o,i] with 

e^/^||Vu5j|L4(n,j,) 0 and Jg \ 0 as e \ 0. 


(60) 


Testing ( [23| with f = Ug — us^ (possible due to (59b)), applying the uniform monotonicity of ^ekPei 
(assumption and ([4^ for p = 4 (compare with the calculation performed in ([4^) gives 


|||e(Mg) - e(M)||^2(Q,^) 

< ?7l|e(M) - e(M5j||^2(f2,j,) + bl|e(Mg) - e{us^)\\\2(^xiT) + ^CudlVug - 

< v\\e{u) - e{usMl^(nT) 

+ / (5elTel(e(Me),Cg,Ze) - 5eWel(e(M4),Ce,Zg)) : (e(Mg) - e(M 5 j)da::dt 

J 

+ e (iVMgpVMg - iVMiJ^VwiJ : (VMg - VM^Jdxdt 

J CIt 
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= r]\\e{u) - e{usj\\l 2 (^n^) 

+ / deWelie{Ue),Cs,Ze) : {e{us) - e{usj) + el^Uel'^Vu^ ■■ {Vue -Vusjdxdt 

J Qt' 


=0 by 


/ deWei(e(usJ,Ce,Zs) I (e(ue) - e{usj) dxdt 

J Q'j’ 

e j : (Vm^ - Vu^J dxdi . 

J Qt- 


(61) 


(*) 


Finally, 


|(*)| < e||Vu5j|i4(n,^)||Vue - VM5j|L4(n,^) 

< (e^^'^WVusAlL^nT)^ (£^^A\"^uA\L^(n^)+e^/A\'^usA\LHnj 


>0 as £\,0 by (|60[> 


<C by Lemma|5.14| 


>0 as £'\0 by l|60[l 


From growth condition (llal, Lemma 5.15 and Lebesgue’s generalized convergence theorem, we 
obtain 

deWei{e{us,),Ce,Ze) deWei{e{u), c, z) in L^{0.t) 
for a subsequence £ \ 0. By ^ u in L°°([0, T]; iL^(n; M")) for a subsequence £ \ 0 (Lemma 


5.15 (hi)) as well as (59a), we also have 

e{ue) - e{us^) ^ 0 in 


as £ \ 0 for a subsequence. Therefore, every term on the right hand side of (611 converges to 0 as 
£ \ 0 for a subsequence. This shows Mg —>■ m in L^([0, T]; iJ^(n; K")) as £ \ 0 for a subsequence 
by Korn’s inequality. 


(ii) Testing (22) with Cg and c and passing to £ \ 0 for a subsequence eve ntual ly shows strong 
convergence Cg —>■ c in L^{[0,T]; (see the argumentation in Lemma 

L s{dtCs)cs dxdt < £\\dtcA\L 2 (nj.)\\ce\\L 2 inj,) -)> 0 as £ \ 0). 


5.10 


and notice that 


(hi) According to Lemma with f = ( = z and /m = (here we choose £m = 1/Af) we find an 
approximation sequence {Cefc} ^ LP{[0,T]]W^’^{Q)) n L°°{Qt) with £fc \ 0 and the properties: 

Ce, ^ z in LP{[0,T];W^’P{n)) as koo, (62a) 

0 < Cefe < ^ck for all k gN. (62b) 

We denote the subsequences also with {ze} and {Ce}, respectively. The desired property Zg 


in LP{[0,T];W^'P{il)) as £ \ 0 follows with the same estimate as in the proof of Lemma 5.11 by 
using the uniform convexity of a; i—>■ \x\p and the integral inequality (24) with C ■= Ce ~ (note 
that {veXe — Ze) = d holds by and (62b)). Indeed, we obtain 


|Vz,-Vzrdxdf 

J 1 ^ 7 " 

< \\dzWe\{e{Ue), Ce, Zg) — a + l3dtzA\L2{[0,T]-,L^{n)) IlCe ~ ^e||L2([0,T];L“(n)) 


bounded 


+ l|V^e|li;L ) live - VzlUp(n^) - / iVz^-^Vz • V(z, - z) dxdt 

V__s/-^ 


bounded 


as £ \ 0 for a subsequence. Her e, we have used z^ ^ z and e “t z in L^([0, T]; L°°(H)) as £ \ 0 
for a subsequence due to Lemma 5.15 and the compact embedding VF^’P(H) ^ L°°(H). ■ 
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Corollary 5.17 The following convergence properties are fulfilled: 


(i) z, ^ z in LP{[0 ,T]-W^’P{Vl)), 
Ze{t) z{t) in a.e. t, 

Zg ^ z a.e. in fix and 
Ze ^ z in W{[0,T]-L'^{il,)), 

(a) Cs ^ c in ([0, T]; H^{Q)), 
Ce{t) —>■ c{t) in a.e. t and 

Cc —> c a.e. in fir, 
as e \ 0 for a subsequence. 


(Hi) Ue ^ u in L^([0,T];_ff^(ri;IR")), 

Ug{t) u{t) in a.e. t and 

Ue ^ u a.e. in fir, 

(iv) /ie ^/i m L^([0,T];7Ji(f])), 

(v) dcWch{ce) dcWchic) in L’^iVlx) 


Now we are well prepared to prove the main result of this work. 


Proof of Theorem 4.6 We can pass to £ \ 0 in (22) and (23) by the already known convergence 
features (see Corollary 5.17) noticing that e\\7Ue\^'^: V^d^dt and e{dtCe)C dxdt converge to 
0 as £ \ 0. We get 


deWe\{e{u), c, z) : e(C) dxdt = 0 


(63) 


for all ( G L^i[0, T]; Wp’‘^(f2; M”)). A density argument shows that (63) also holds for all ( G A^([0, R”)). 

Writing (21) in the form 

/ (cg — c°)9tC dxdt = / dxdt, 

J Ot-- J Ot-- 

by only allowing test-functions C G ^^([0, T]; 77^(11)) with dt( G and ({T) = 0, we can also pass 


to £ \ 0 by using Corollary 5.17 


To obtain a limit equation in (24) and ( |25[ ) , observe that 

[5^Wei(e(u£),Ce,Ze)]+ ^ [5^Wei(e(u),c,z)]+ in L^iVtr), 

X{z,=o}^X, inL^inr) 

for a subsequence £ \ 0 and an element y G L°°(f2T)- Setting r := —x[i92lTei(e(u), c, z)]+ and keeping 
(27) into account, we find for all C G L°°{Qt)' 


dxdt 


r( dxdt 


(64) 


for a subsequence £ \ 0. Thus, we can also pass to £ \ 0 for a subsequence in (24) by using Lebesgue’s 
generalized convergence theorem, growth condition (A5), Corollary 5.17 and (64). Let ^ G L°°([0,T]) 
with ^ > 0 a.e. on [0,r] be a further test-function. Then, (25) and (27) imply 

0> / ( f re(t)(C- 0 e(t))dxV(t)dt= [ r sic - Zs)£, dxdt 

Jo \Jn J Jqt 

[ r{C — z)Cdxdt= f f f r(t)(C — ^(t)) dx^ ^(t) dt. 

J Jo \Jq, j 

This shows r(t)(C — z{t)) dx < 0 for a.e. t G [0,T]. 

It remains to show that (26) also yields to a limit inequality. First observe that (26) implies: 

£eiqsit 2 ))+ [ aiZeih) - Zeit 2 )) dx + f / /3|9t^eP + | V/ie^ dxdt - 4(ge(ti)) 

J Q J ti J Q 
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< 


/7 

/ ti 7 r2 


deWe\{e(ue),Ce,z^) ■. e{dth)(ixdt + e f f \Vue\^Vue : Vdtbdxdt. 

J t\ «/ O 


(65) 


To proceed, we need to prove e |Vue(t)|‘^dx —>■ 0 as e \ 0 for a.e. t € [0,T]. Indeed, testing (23) with 

^ := Ug — b gives 

£ / |Vue|^dxdt = e / iVuepVue : V6da;dt — / deWe\{e{ue),Ce, Zg) : e{ue — b) dxdt. 

J J J 

We immediately see that the first term converges to 0 as e \ 0. The second term also converges to 
0 because of deWei(e(u),c, z) : e(u — b)dxdt — 0 (equation (63)). This, together with Corollary 


5.17 proves £e{qe{t)) —>• £(q{t)) for a.e. t e [0,T]. In conclusion, we can pass to e \ 0 in (65) for a.e. 


0 < fi < ^2 < T by Corollary |5.17| together with Lebesgue’s generalized convergence theorem, growth 
condition (|A2|), (11a) and (A6) as well as by a sequentially weakly lower semi-continuity argument for 


/3|9(Zep dx and for iVTiTf dx. 
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